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We develop a theory of non-linear cosmological perturbations at superhorizon scales for a scalar 
field with a Lagrangian of the form P{X,(j)), where X — —d'^cpdfi.cj) and cj> is the scalar field. We 
employ the ADM formalism and the spatial gradient expansion approach to obtain general solutions 
valid up to the second order in the gradient expansion. This formulation can be applied to, for 
example, DBI inflation models to investigate superhorizon evolution of non-Gaussianities. With 
slight modification, we also obtain general solutions valid up to the same order for a perfect fluid 
with a general equation of state P = P{p)- 
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I. INTRODUCTION 



Generation of primordial fluctuations during inflation is one of the most interesting predictions of quantum field 
theory. Indeed, those quantum fluctuations are considered as seeds of the large scale structure of the present universe, 
and this picture has been accepted by many researchers as a standard scenario. The cosmic microwave background 
(CMB) temperature anisotropy found by COBE [l| was perfectly consistent with the predictions of the linear theory of 
^f,- cosmological perturbations Q- In particular, the primordial fluctuations are nearly scale invariant and consistent 
_j' with Gaussian statistics 

The recent more accurate observation by WMAP, however, has revealed deviation from exact scale invariance, 
with a slight red tilt [H. Moreover, there is a good possibility that deviation from Gaussianity can be detected 
by the future experiments such as PLANCK 0]. With those current and future precision observations, deviation 
from the exact scale invariance and Gaussianity can be a powerful tool to discriminate many possible inflationary 
models. Especially, non-Gaussianity of primordial fluctuations has recently been a focus of much attention by many 

authors [1, i M M M M JA, la, M E Ha, la, 20, 21, 22. 2^.m.2' 



In order to parameterize the amount of non-Gaussianity of primordial perturbations, commonly used is the non- 
linear parameter Jnl- This is related to the bi-spectrum of the curvature perturbation, and is defined as [T6| 



00 , o 

P.- C(x)-Cg(x)--W^(x), (1.1) 



where Cg is the curvature perturbation on a uniform density hypersurface and satisfies linear Gaussian statistics. On 
the observational side, current bounds on the parameter /jvl by WMAP five years @ are —9 < Jnl < 111 for the 
■ local form of bi-spectrum and —151 < Jnl < 253 for the equilateral form. B y fu ture experiments such as PLANCK 
01, it is expected that non-Gaussianity of the level I/atlI^S can be detected |16| . 

On the theoretical side, although the non-Gaussianity from the standard single-field slow-roll inflation is suppressed 
by slow-roll parameters [l^ and is too small to be detected in near- future experiments, many new types of 
inflationary models predicting large non-Gaussianity {/nl greater than unity) have been proposed. There are at 
least two known mechanisms to generate large non-Gaussianity: isocurvature perturbations and non-canonical kinetic 
terms. Intriguingly, bi-spectrums for these two mechanism typically have different shapes: the so called local type 
for non-Gaussianity from isocurvature perturbations and the equilateral type for non-Gaussianity from non-canonical 
kinetic terms. 

As for the former case, a typical example is the curvaton scenario [1^, [13, [H, H^, where a light scalar field, 
called curvaton, is responsible for isocurvature perturbations during inflation. Another example is multi-field inflation 
models, where generation of isocurvature perturbations is generally expected because of the existence of more than one 
fields in the inflaton sector In these examples, isocurvature perturbations generated during inflation 
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are later converted to curvature perturbations, and this process is important for large non-Gaussianity. Now let us 
consider the later case, where non-canonical kinetic terms are responsible for l arg e non-Gaussianity. Examples of this 
type include k- inflation [13, [IH, ghost inflation [Hjls^l and DBI inflation [sljIlHl. In fact, Weinberg [H] has recently 
shown that the leading corrections to the Gaussian correlations of curvature perturbation are solely of the k-inflation 
type (including DBI inflation), except for ghost inflation. In k-inflation and DBI inflation, large non-Gaussianity is 
expected when the non- linear nature of the non-canonical kinetic action becomes significant. This happens when the 
sound speed of perturbations is sufficiently smaller than unity [STf . 

On the other hand, in ghost infiation large non-Gaussianity, such as I/nlI — 80, is always expected unless fine- 
tuned. The reason is that non-linear terms in the low energy effective Lagrangian is suppressed only by fractional (not 
integer) powers of energy scales. To be more precise, the leading non-linear term is suppressed only by (ff/M)^/^, 
where H is the Hubble expansion rate during inflation and M is the cutoff scale of the low energy effective theory. 
This is confirmed by a simple scaling analysis and, thus, should be robust. After using the COBE normalization 
Sp/p ~ 10"^, which implies H/M ~ 10~*, this fact leads to the prediction of large non-Gaussianity [s^ if the 
dimensionless coupling constant for the leading nonlinear term is set to be order unity, i.e. if we do not fine-tune the 
theory. 

For these reasons, non-Gaussianity is one of the most powerful tools to distinguish models of inflation with combi- 
nation of the future observations. Thus, to quantify the non-Gaussianity and clarify its observational signature, it is 
important to develop a theory that can deal with nonlinear cosmological per turbations. There are couple of methods 
to tackle this problem. One is a second-order perturbation theory l38j . Another is based on spatial gradient 
expansion (isl. [soj. While the former mainly deals with primordial perturbations up to around the horizon exit, the 
later deals with classical evolution after horizon exist. Thus it is important to develop both methods and to use them 
complcmentarily. 

Closely related to the gradient expansion method, cosmological perturbations on supcrhorizon scales have been 
studied extensively in the so-called separate universe approach or 5N formalism li^ . [4l[ . Actually, these approaches 
are essentially the leading order approximation to the gradient expansion p^ . Il3l |. Including these, many of the 
previous studies were confined to the leading order approximation to the gradient expansion. However, higher order 
corrections to the leading order results can be important to get more detailed information about non-Gaussianity. 
One good example is the case studied by Leach et al [42| . They considered linear perturbations in single-field infiation 
models and supposed that there is a stage at which slow-roll conditions arc violated. It has been then shown that, 
due to the decaying mode, the O(e^) corrections in spatial derivative expansion do affect the evolution of curvature 
perturbations on superhorizon scales [i^. A similar situation for DBI infiation was considered by Jain et al [i^. In 
these situations, it is expected that non-Gaussianities should also be affected by the order O(e^) effects. However, the 
linear perturbation theory is not capable for calculation of non-Gaussianity. Thus, it is necessary to develop nonlinear 
theory of cosmological perturbations valid up to O(e^) in the spatial gradient ex pan sion. 

Gradient expansion formalism has been developed and used by many authors [l3l . [2^ . [25l . [sol . [iol . [4]| . Formulation 
valid up to O(e^) was developed, for example, by Tanaka and Sasaki for a universe dominated by a perfect fluid with 
a specific equation of state P/p = const j24| and that dominated by a canonical scalar field (25[. However, as far 
as the authors know, those works have not extended to a perfect fluid with general equation of state P = P{p) nor 
to a scalar field with non-canonical kinetic action, which is essential for the second type of mechanism of generating 
non-Gaussianity. 

The purpose of this paper is to fill this gap. Namely, we shall develop a theory of nonlinear superhorizon perturba- 
tions valid up to the order O(e^) for a scalar field with non-canonical kinetic action and a perfect fiuid with general 
equation of state. 

This paper is organized as follows. In Sec. [Ill we introduce a non-canonical scalar field as our model and express it 
in a perfect fluid form. In Sec. IIIII we shall develop a theory of nonlinear cosmological perturbations on superhorizon 
scales and explain our formulation in details. The following Sec. IIVI is devoted to some details of obtaining a general 
solution. We then study some specific examples in Sec. [V] Sec. IVII is devoted to a summary of this paper and 
discussion. In Appendix, we give our result in a perfect fiuid system. 

II. SCALAR FIELD IN A PERFECT FLUID FORM 



Throughout this paper we consider a minimally-coupled scalar field described by an action of the form 

/ = / d^Xy^P{X,(P), 



(2.1) 
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where X = ~g^'^df^(j)d^(l>, and suppose that —g^'^di.cj) is timehke and future-directed. The equation of motion for (f> is 

^ d^{V^Pxd''<l>)+P^ = Q, (2.2) 



where the subscripts X and </> represent derivative with respect to X and (j), respectively. The stress energy tensor of 
the scalar field is shown to be a perfect fluid form: 

T^, = 2Pxdf,^d,(t> + Pg^, = (p + P)u^u, + Pg^„ (2.3) 

where 

p{X,<l,)^2PxX-P, M^ = -M. (2.4) 

V A 

Note that u^u^ = —1. As far as d^^cf) ^ 0, the equation of motion (|2.2p is equivalent to the conservation equation 
The following relation among first-order variations of P, p and (j) will be useful in the analysis below. 

5P = cl5p + pT5(j), (2.5) 

where 

Cs = ^p . r = i (P^ - c^p^) . (2.6) 

Note that is the speed of sound for the gauge invariant scalar perturbation in the linear theory [lO] • 

III. FORMALISM 

In this section we shall develop a theory of nonlinear cosmological perturbations on superhorizon scales. For this 
purpose we employ the ADM formalism and the gradient expansion in the uniform Hubble slicing. 

A. ADM decomposition 

In the (3 + l)-decomposition, the metric is expressed as 

ds^ = gf.udx'^dx'' = ^a^dt^ + 7y(dx* -I- P'dt){dx^ + (3^dt), (3.1) 

where a is the lapse function, /3* is the shift vector and Latin indices run over 1,2,3. Since a and represent 
gauge degrees of freedom for diffeomorphism and appear as Lagrange multipliers in the action, the corresponding 
equations of motion leads to constraint equations. Contrary to a and /3, components of the spatial metric jij 
are dynamical variables (subject to the constraint equations) and the corresponding equations of motion are called 
dynamical equations. In what follows we shall express the dynamical equations as a set of first-order differential 
equations with respect to the time t. For this purpose we introduce the extrinsic curvature Kij defined by 

K^j = - ^ {^a^, - D.P, - D.Pi) , (3.2) 

where D is the covariant derivative compatible with the spatial metric ^ij . For the stress-energy tensor in the perfect 
fluid form (|2.3p . we define the 3-vector ii* as = /u^. Hereafter, we shall use 7^ and its inverse 7'^ to raise and 
lower indices of D, v, (3. Then we can express and in terms of a, 13"^ and v^: 



[a^^ («,+/?fc)(^'' +/?')] 



-1/2 



= u"(w,+ft)- (3.3) 
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The conservation equation VmT'^ = is 



1 



--dt 



\/det7(p + P)wu^ + [(p + P)^?;^'?/^] -ad^P - {p + P)w [wd^a - UjD.P^] 



where w = au^ . All independent components of the stress energy tensor are conveniently expressed as 

E = Tf^^n^n" = {p + P)w^ - P, 
Ji = -Tf^iu'^ = {p + P)wui, 
Sij = Tij = {p + P)uiUj + P"iij, 

where is the unit vector normal to the constant t surfaces and is given by 



(3.4) 



(3.5) 



(3.6) 



The Hamiltonian constraint, corresponding to the equation of motion for a, is 

R + K^ - KijK'^ = 2k^E, 



(3.7) 



where — SttGn, R = R[j] is the Ricci scalar of the spatial metric 7 and K = K^^. The momentum constraint, 
corresponding to the equation of motion for /?% is 



The dynamical equations, decomposed into the trace part and the tracelcss part, are 



(3.8) 



- W^K,, + KK, 



D^DjU 7.y 



R - d^K + A'2 



(3.9) 



where = n^d^, V_l = n^V^, = -f^^ DiDj and Rij = Rij['j] is the Ricci tensor of the spatial metric 7. 

In addition to the standard ADM decomposition briefly reviewed above, we further decompose the spatial metric 
and the extrinsic curvature as 



7y = aV7ij-, 



(3.10) 



where a(t) is the scale factor of a fiducial Friedmann background (specified later) and the determinant of 'jij is 
constrained to be unity: det7y = 1. The first-order equations for the spatial metric jij) are deduced from the 
definition of the extrinsic curvature (13.21) as 



d±'>P _^ dta 
tp 2aa 



-K 



a 



a 



(3.11) 
(3.12) 



The first-order equations for the extrinsic curvature {K, Aij) arc obtained from the dynamical equations p.9p as 



.... 

3 



{S\+E), 



R 



1 



5,, 



(3.13) 



(3.14) 
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where A.^^ = j^'^Aij and A^'^ = '^'^^ A^^ . The Hamiltonian and momentum constraints are, respectively, 

i?+ -i'^ly = 2k^E, (3.15) 

~'^d,K + D,Ar> = k'J,. (3.16) 

The conservation equation is 

{dt + v'di) p + {dt [{ai^'^fw] + [{axl,''fwv'] } = 0, (3.17) 

^ dt [{ailj'^fip + P)wu.,] + Dj [{p + P)wv^u,] + d,P + (p + P){w^d^a - wujD,J3^) = 0. (3.18) 



(07/^2)3 

Throughout this paper we adopt the uniform Hubble slicing 



K=-iH{t), H{t) = —. (3.19) 

a 



Substituted into p.ll|) . this imphes that 



B. Gradient expansion: basic assumptions and order estimates 

In the gradient expansion approach we introduce a flat FRW universe (a(t), (l)o{t)) as a background and suppose 
that the characteristic length scale L of perturbations is longer than the Hubble length scale 1/H of the background, 
i.e. HL ^ 1. Therefore, we consider e = 1/{HL) as a small parameter and systematically expand our equations by 
e, considering a spatial derivative acted on perturbations is of order 0(e). 

The background flat FRW universe {a(t), (t>o(t)) satisfies the Friedmann equation and the equation of motion 

= ^Po, ^dt {a^Poxdt(l)o) ~ Po^ = 0, (3.21) 

where po = p{Xo, 0o), Pqx = Px{Xo, 0o), -Po0 = -P^l^o, M, and Xq = {dt4>of ■ 

Since the FRW background is recovered in the limit e — > 0, we naturally have the estimates 

v' = 0{e), /3* = 0(e), (3.22) 

and dt^ij = 0{e). Actually, following the arguments in rcfs. (2^. [25|. we assume a stronger condition 

dt% - O(e'). (3.23) 

This assumption significantly simplifies our analysis and, we believe, still allows many useful applications of the 
formalism. On the other hand, we consider -0 and (without derivatives acted on them) as quantities of order 0(1). 

We can estimate orders of magnitude of various quantities by using the above assumption and the basic equations. 
First, ([XT^ imphes that 

A,,=0{e^). (3.24) 

Substituting (|3?24l) into (|3T6l) we obtain Ji = 0{e'), or 

z;, + A = 0(e3). (3.25) 

For the scalar field system, this is expressed as dii: — 0{e^), where tt = <p — (j)o. By absorbing a homogeneous part of 
TT into 00 (and modifying a{t) accordingly), we obtain 



TT = 0{e^). 



(3.26) 
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Combining (|3.25p with the first equation in (|3.3p . we obtain = l/a + O(e^), or 

w (=au°) = l + 0(e6). (3.27) 
This impUes that E ^ p + O(e^). Thus, from ([STTS]) we obtain 



and (|2.5|) imphes that 



FinaUy, (|3.17p imphes that 



and thus wc obtain 



S ( = ^^) =0(e2), (3.28) 
Po 



p (= P - Po) = 0{e^). (3.29) 
dt^j = 0{e^), (3.30) 
X = 0(e2) (3.31) 



from the uniform Hubble shcing condition ()3.20p . 

In summary, we have the following estimates (including assumptions): 

V^ = 0(1), 7y=0(l), v'=0{e), /3'=0(e), 

X = 0(e^), A,,-0(e2), 5 = ©(e^), tt = 0(6^), p = 0(e2), 

9*7., =0(e2), 9t^ = 0(e2), t;' + = ©(e^), w = l + 0{e''). (3.32) 

C. Leading order equations 

Substituting the order of magnitude shown in (|3.32|) into the conservation equations p.lTp and (|3.18p . we find 

PodtS + (po + Po) (^6^ + Dy^ + 3H{p - PoS) = 0{e^), (3.33) 
\dt [a3(po + Po)u..] + d, [p + {po + Po)x] = O(e^). (3.34) 
The Hamiltonian and momentum constraint equations give 



R-2K^a^iP^P(^5 + 0{e^), (3.35) 



(^^i/) = K\p^ + P^)i;^u, + 0{e''), (3.36) 



where R = R['f\ is the Ricci scalar of the normalized spatial metric 7y, D is the covariant derivative compatible with 
7ij, = DiDj, and 7'-' is the inverse matrix of 7^^. The evolution equations for the spatial metric give 

^ 3i/x + A:/3'+0(e4), (3.37) 

V' 

{dt - P''dk)% = -2A^j + 7,fc9,/3^- + %k^^^3'' - ^7y5fe/3'= + 0{e^), (3.38) 
while the evolution equations for the extrinsic curvature give 

dtA,j + ?.HA,j = {r,j - \Rl.?j + O(e^), (3.39) 

P+iPo + Po)x = -lpo5 + 0ie^), (3.40) 
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By using (|3.37p . (|3.40p and the background conservation equation dtpo + 3H{po + Pq) = 0, a single equation for 5 
is easily obtained from p.33p . 



Using (|3.40p again, p.34p is simplified to 



(3.41) 



(3.42) 



It is intriguing to note that we have not yet specified the form of p. Therefore, eqs. p.33p - p.42p hold as far as 
the stress-energy tensor is of the perfect fluid form and p = O(e^), regardless of whether the stress-energy tensor is 
provided by a scalar field, radiation, dust, or any other sources. 

The form of p for the scalar field system is specified by the relation ()2.5p as 



P = po{cl„5 + To7r)+Oie^), 



(3.43) 



where c^q ~ Pqx / i'^Pox x Xq + Pqx) and To = (Po^ ~ '^mP04>)/ Po- We can obtain another equation relating p and tt, 
by expanding p as p = PoxiX — Xq) + Po,pT^ + O(e^), where X — Xq = 2[dt4>{)dt'n — xXq) + O(e^). Actually, this 
equation can be interpreted as a first-order equation for tt. Using (|3.40p . pq + Pq = 2PqxXq and the background 
equation of motion (j3.2ip . we can rewrite this equation for tt as 



(po + Po)a^ 



5, 



-Pod + Oie* 



(3.44) 



IV. GENERAL SOLUTION 

Having written down all relevant equations up to the order O(e^) in the gradient expansion, we now seek a general 
solution. 



A. Leading order solutions 



First, V = 0(1) and dtip = 0{e^) imply that 



V' = i(°)(.T'=) + 0(e2), 



(4.1) 



where L'^'^\x^) is an arbitrary function of the spatial coordinates {x''} (k ~ 1,2,3). Hereafter, the superscript (n) 
indicates that the corresponding quantity is of order 0(e"). Similarly, jij = 0(1) and dtjij = 0{e^) imply that 



(4.2) 



where fl^\x^) is a (3 x 3)-matrix with unit determinant whose components depend only on the spatial coordinates. 



With these expressions, the right hand side of (|3.39p is calculated as 



(4.3) 



where 



{my 



2 ( 2d, In InL(o) - D^°^ d':!'^ In^ 



-|/« (2a,lni(°)5,lnL(°) -^f ^(°)lnL(")) fjf 



(4.4) 



ffi^ is the inverse matrix of f^^\ = and i?'"' ~ are Ricci tensor and Ricci scalar of the Oth-order 



spatial metric f^j , and I?'") is the covariant derivative compatible with f^j . Note that f(^-^F^^ — by definition. 
Thus, we obtain 



A, 



'it) 



F!f\x') I a{t')dt' + Cl^\x'^) 



(4.5) 
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where C^^(a;'=) is a symmetric matrix whose components depend only on the spatial coordinates and which satisfies 
Next, the equation (|3.4ip for 5 is easily solved as 



5 = 



PoitWit) 



(4.6) 



where p^a^ is a constant and Q^^''{x'^) is an arbitrary function of the spatial coordinates. With this expression for S, 
(IPS]) gives 



Ui 



p*a^ 



Z[po{t)+Po{t)]a^{t) 



d^Q^^\x'') / a(Odi' + cf' (x'=) 



(4.7) 



where Cp''(a;'"') is an arbitrary function of the spatial coordinates. 

The 'constants' of integration L(°)(x''), f^\x''), c\f(x^), Q^'^\x'') and cf'\x^) are not independent. Indeed, by 

solving the Hamiltonian and momentum constraints (|3.35p and (|3.36p . Q^^) and C^^^ are expressed in terms of other 
integration 'constants' as 



1 



3 



R 



{my K^p^ai 



pk f,(0) 



+ 0(e^). 



(4.8) 



Until now, we have not used either p.43p or (|3.44p . Therefore, the general solutions presented above are valid not 
only for the scalar field system but also for radiation, dust or any other sources, provided that the stress-energy tensor 
is of the perfect fluid form and that p = 0{f^). 

We now use p.43p and (|3.44[) to proceed further. It is easy to integrate (|3.44[) to give 



P*aldt<j)o 



Q(2)(2;fe) / a(f')rf^' + n(2)(x^) 



to 



3(po + Po)a^ 

where n(2)(a;'^) is an arbitrary function of the spatial coordinates. By using (|3.43p and (|3.40p we obtain 



(4.9) 



Hpo + Po)a^ 



1 + 3c% 



Po^odt(l)o 
{po + Po)a Jt„ 



a{t')dt' 



(po + Poja 



(4.10) 



Note that, since Ui = ~di(j)/VlC = —diT:/dt<po + 0{e'), the integration 'constant' Ii^'^\x'') is related to cf\x^) as 



(4.11) 



Thus, Cj for the scalar field system does not include a transverse part. 



B. Solution up to O(e^) 



Solutions obtained so far are correct up to leading order in the gradient expansion. Among them, the spatial metric 
■0 and have been obtained only up to 0{\) while all other variables are correct at least up to O(e^). In this 
subsection we seek 0(6^) corrections to V' and 7ij. For this purpose it is convenient to specify the shift vector /3* more 
accurately than indicated by (j3.32p : in this subsection we set 

f3 = 0{e-). (4.12) 

With this gauge choice, (|3.38p is reduced to 



^a^J=-2A^,+0{€^), 



(4.13) 
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and thus results in 



dt' 



<^'it') 



t' rt 
a{t")dt" + Cl-\x'') 



to 



(4.14) 



where we have absorbed a new integration 'constant' into the 0-th order integration 'constant' fl'^\x^). Again, we 
have not used either p.43p or (|3.44p to derive this result. Therefore, the general solution (|4.14p is valid not only for 
the scalar field system but also for any other sources, provided that the stress-energy tensor is of the perfect fluid 
form, that p = O(e^) and that the additional gauge condition (|4.12p is imposed. 

In order to seek the O(e^) correction to i;, we specialize to the scalar field system since we will need the solution 
(|4.10p for X, which was obtained by using (|3.43p . The evolution equation p.37p with (|4.ip leads to 



In 



[ ^ 1 




_L(0)(a;'=)_ 


2 A 



H{f)x{t')dt' + 0{e% 



(4.15) 



where we have absorbed a new integration 'constant' into the 0-th order integration 'constant' L^^\x'^). We can 
substitute (|4.10p to the right hand side of this equation to complete the procedure. 



C. Summary of the solution 

In summary we have obtained the following solutions in the gradient expansion for the scalar field system. 



1 



1 



R 



dt(f>a 



1 



d, ( R 
R 



(i(0))4j(0)l / a(i')rft' + C(2) +0(e5), 



6K2(po + Po)a2 

+ 0(6^), 



(i(0))4_^(0) / a{t')dt' + C^^^] +0{e^), 

J J to 

1 , q 2 POrO^t^O f^,^^^, 



{po + Po)a Jto 



a{t')dt']R (l(°))V^°' 



poTQdt(t)o ^(2) 
(po + Po)a 



(0) ,1^(2) l' dt' f* ^(^„)^^„^^(2) dt' 



7., = n'-2 



A, 



toO^'{t')Jto ' ' Jtoa^t') 



^ ' f}?^ I a{t')dt' + C\f ) +0(6^), 



(4.16) 



where C^^) this expression is related to II^^) in (|4.10p as C'^) — 2k'^ p^a1Xl^'^\ and F>j ' is defined by (|4.4p . and 

'constants' of integration L^'^\ f^^\ C^^) and C|J' depend only on the spatial coordinates {x^} {k = 1,2,3), and 
satisfy 



? = ff. det(if)) = l, 

(2) _ ^(2) „y ^(2) _ „ 



a 



(4.17) 



Here, /^'^^ is the inverse matrix of /^'j'' and Z)^^) is the covariant derivative compatible with f^^- 
Note that the gauge condition (|4.12p is unchanged under purely spatial coordinate transformations 

x' ^ x" = P{x^). 



(4.18) 
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Thus, the 0-th order spatial metric fj^j^ includes 3 gauge degrees of freedom. Therefore, the number of degrees of 
freedom included in each 'constant' of integration is 



i(0) 

AO) 
Jij 



1 scalar growing mode = 1 component , 

2 tensor growing modes = 5 components — 3 gauge , 

1 scalar decaying mode — 1 component , 

2 tensor decaying modes = 5 components — 3 constraints . (4.19) 

V. EXAMPLES 

Until now. we have not specified a form of the function P(X,(j)). In this section, we shall consider some specific 
examples. 

A. Scalar with shift symmetry 

Let us consider the case where the Lagrangian P depends on the scalar field only through X, i.e. P ~ P{X). This 
case corresponds to a scalar with shift symmetry and is often considered in models of k- inflation [3l| . 
The sound speed c^g ^'^'^ ^o in this case are 

4- — , ro = o, (5.1) 

Pox 

and the sound speed agrees with the adiabatic sound speed defined by SP — c^Sp. From this fact, it is expected that 
dynamics of the scalar field system should somehow resemble that of a perfect fiuid with equation of state P = P{p)- 
This expectation turns out to be essentially correct. In fact, in Appendix we see that these two systems have essentially 
the same general solutions up to O(e^) in the gradient expansion except for the following one difference. The 3- velocity 
Ui for the scalar field system does not include transverse mode but that for the perfect fiuid in general does. The 
evolution of other quantities such as i5, y, jij and Aij are the same. 

Note that ghost condensation 0, [4^, also has a similar low energy effective Lagrangian but includes terms like 
(V^7r)2 as well, where V is the spatial gradient and tt is perturbation of the scalar field. Indeed, those additional terms 
play important roles in generation of primordial density perturbations [s^ and infra-red modification of gravity [4^ . 
In other words, as easily seen by doing a proper scaling analysis, a Lagrangian of the form P{X) without those 
additional terms can not describe ghost condensation in general even at low energies. However, if we are interested in 
a situation where the the Hubble scale during infiation is longer than the scale of IR modification of gravity then we 
can safely use the present formalism to investigate classical evolution of cosmological perturbation at superhorizon 
scales. However, if the scale of IR modification is longer then the present formulation is not valid. In this case we 
probably need to extend the present formulation to include the higher derivative term. More detailed investigation 
will be considered in future publication. 

B. Canonical scalar 

Next let us consider a canonical scalar field, i.e. 

P(X,</>) = ^-y(</>), (5.2) 
where V{(f>) is a potential. In this case, we can obtain 

4 = 1, r„ = -^, (5.3) 

Po 

where V^o = dV {(l)o) / d(l)o . The general solution can be read off from (|4.16p . In particular, x is given by 



X 



'303a3 



2 r / /■* 



to 



2a0o + V^o / a{t')dt' Q^^) ^^k^ ^ y^^^{2) ^^k^ 



+ Oie-), (5.4) 



where (po = dtcjjQ. If we set Q'^'^\x'^)p^a1 ~ — 3(/)*C(.t''') and \l'^'^\x^)p^a1 ~ — 30,£'(a;'^), then this reduces to eq. (C29) 

of m. 
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C. DBI scalar 



We now consider a scalar field described by the so called DBI action. This kind of scalar is considered as an inflaton 
in an interesting class of inflationary models called DBI inflation [sj, [11] . In a simple case the inflaton scalar field 
corresponds to the radial position of a D3-brane in a warped compactification. 

For a warped throat background, the 10-dimensional metric takes the following form 

ds^ = h^{p)7jf,,dx^'dx'' + h-^{p) (dp' + p''gi^ldx"'dx''^ , (5.5) 

where h{p) is a warp factor, x'^ {p = 0, 1, 2, 3) are external 4-dimensional coordinates, p is the radial coordinate in the 
warped throat, and (to = 5,6,7,8,9) are the internal 5-dinicnsional angular coordinates. Considering the radial 
position p of a D3-brane in this background as a scalar field in the external 4-dimensional spacetime, its dynamics is 
described by a Dirac-Born-Infcld (DBI) action plus a Chern-Simons term and additional potential terms. Thus, the 
Lagrangian P{X, (p) is 

P{X, (f) = -T(0) Vl - X/T{^) + T(0) - V{4>) , (5.6) 

where </> = T^^^p is the scalar field describing the D3-brane position, T(0) = T'ih'^{(j}) is the warped brane tension, 
and V{4)) is the inflaton potential. Here, T3 is the D3-brane tension. The Lagrangian (|5.6p also applies to the case 
where the inflaton cf) is the radial position of a wrapped D5- or D7-brane (47l |. 
For the Lagrangian (|5.6p we obtain 



1 1 

ro = — 



7 V 7^ 



(5.7) 



where T^q = dT {(f)Q) / dcfiQ and 7 = 1/ y^l — Xo/T((/)o). In a non-relativistic limit (7 ~ 1), c^q and Fg reduce to those 
for the canonical scalar field discussed in subsection IV Bl 



VI. SUMMARY AND DISCUSSION 



We have developed a theory of nonlinear cosmological perturbations on superhorizon scales for a scalar field de- 
scribed by a Lagrangian of the form P{X, (/>), where X — —3^09^0 and </> is the scalar field, and also for a perfect fluid 
with a general equation of state P = P{p). The general solutions valid up to the order O(e^) in the spatial gradient 
expansion have been presented in subsection IIV CI for the scalar field system and in Appendix \^ for the perfect fluid. 

This formalism can be applied to many interesting circumstances. Some particular examples have been listed in 
Sec. |Vl including a scalar with shift symmetry, a canonical scalar and a DBI scalar. (As shown in Appendix [A] it 
can be applied also for a perfect fluid with general equation of state P = P{p) ) Thus, the formalism can be used 
to investigate superhorizon evolution of nonlinear cosmological perturbations in k-inflation and DBI inflation. As 
explained in Sec. HI non-Gaussianity can be affected by order O(e^) corrections if, e.g. there is a stage which violates 
the slow roll conditions. 

For ghost inflation, applicability of the present formulation seems a bit subtle as briefly discussed in subsection IV Al 
If the Hubble scale is longer than the scale of IR modiflcation of gravity then we can use the present formalism to 
investigate superhorizon perturbations. On the other hand, if one is interested in a situation where the scale of IR 
modification is longer then the present formulation is not valid. We hope to address this issue in more detail in future 
publication. 

Recently, models of multi-field DBI have been also studied to investigate large non-Gaussianity by [13, [2^ . Here 
our formulation has been developed in a single scalar field, however, we also plan to extend it to a system of multi-field 
scalar in the future. 
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APPENDIX A: PERFECT FLUID WITH P = P{p) 



111 this appendix, we shall consider a universe dominated by a perfect fluid with a general equation of state P — P{p)- 
The stress-energy tensor is given by 

T^, = {p + P)u^u, + Pg^, , P = P{p) . (Al) 

The formulation in this case is similar to that presented in Sec. IIIII and Sec. IIVI for a scalar field. Main differences 
come from the form of the pressure perturbation p. 

As in Sec. IIIII and Sec. IIV( we employ the ADM formalism and the gradient expansion in the uniform Hubble 
slicing. In the ADM formalism, the dynamical equations are given by (|3.9p . with the two constraint equations (|3.7p 
and (|3.8p . We further decompose the spatial metric and the extrinsic curvature as (|3.10p . Then we obtain the first- 
order equations (|3.11[) and (|3.12p for the spatial metric (■0,7^), and (|3.13p and (|3.14p for the extrinsic curvature 
{K,Aij). The Hamiltonian and momentum constraint are, respectively, (|3.15p and p.l6p . The conservation equation 
V^T/f = is given by ([STT]) and |37T8)) . We adopt the uniform Hubble shcing as ((3?T9)) . 

Next, we will employ the gradient expansion. In this approach we introduce a flat FRW universe (a(i), po(i)) as a 
background and suppose that the characteristic length scale L of perturbations is longer than the Hubble length scale 
\/H oi the background, i.e. HL ^ 1. Therefore, we consider e = 1/{HL) as a small parameter and systematically 
expand our equations by e, considering a spatial derivative acted on perturbations is of order 0{e). The background 
flat FRW universe {a{t), po{t)) satisfies the Friedmann equation and the conservation equation dtPo + 3H{po + Pq) = 0, 
where Pq = P{pq). We can estimate order of magnitude of various quantities by using the assumption (|3.22p and 
(|3.23p . As a result, we have the estimates shown in p.32p except that the condition tt = O(e^) does not exist in the 
present case. Then substituting the order of magnitude (|3.32p into the basic equations gives leading order equations 
as (|333ll - ([3:40ll . 

By using these equations and background conservation equation, equations for i5 and Ui are easily obtained as p.4ip 
and (|3.42p . They are easily solved as (|4.6p and (|4.7p . The traceless part of the extrinsic curvature Aij is solved by 
using the leading part of ip and jij, ()4.ip and ()4.2p . as (j4.5p . The 'constants' of integration are not independent but 
arc related to each other by the two constraint equations as (|4.8p . However, in the present case there is no relation 

like (j4.1ip . and then C^'' contains transverse component as well as the longitudinal component. 



For P = P{p), perturbation of the pressure p = P — Pq is given by 



dPo 
dpo ■ 



By using (|3.40p . we obtain the solution for x as 

(l + 3ego)i?[(L(°))VW] 
^ 6K^po + Po)a^ 



(A2) 



(A3) 



In order to obtain a general solution valid up to O(e^), we further seek O(e^) corrections to t/j and 7^. For this 
purpose we adopt the gauge condition (|4.12p . With this gauge choice, the 0{e'^) part of 7ij can be obtained by using 
the solution of Aij as (|4.14p . Similarly, by using the solution of x (|A3p . the O(e^) correction to ip is obtained as 
(I4.15p . provided that x in the present case is given by (|A3p . 

In summary we have obtained a general solution valid up to O(e^) for the perfect fluid with P = Pip). 



S = 



X 
7y 



2K?pao? 



R 



&n'-{po + Po)a^ 
1 + 3c2„ 



d,R 



(i(0))4j(0) / a{t')dt' + Cr ]+0{e>), 



to 



6K^{po + Po)a^' 
L(0) [l + \j' H{t')x{t')dt^+0{e% 



cm 



2 F„ 



(2) 



dt' 



I- 0?{t') 



to 



^{F}Pla(t')dt' + C^' 



0[e% 



* a{t")dt" + ^* 



(2) 



(A4) 
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where 



C, 



(3) 



pk f,(0) 



(2) 



+ Oie% 



(A5) 



Here F^p is defined by 
{x''} [k = 1, 2, 3) and satisfy 



and 'constants' of integration L'^"^ and depend only on the spatial coordinates 



G 



p(0) 
(2) 



f(0) 

ji 1 

(2) 



dct(/^') 

(2) 



= 1, 

0. 



(A6) 



Here, J^^L is the inverse matrix of jf^ and _D(°) is the covariant derivative compatible with ^fl'' ■ Compared with the 



(0) 



'(0) 



solution (|4.16p for the scalar field system, only differences are: (i) tt does not exist; (ii) To = 0; and (iii) cf^ includes 
not only a longitudinal component but also transverse components. 

The number of degrees of freedom included in each 'constant' of integration is 



c. 



1 growing adiabatic mode of density perturbation = 1 component 

2 tensor growing modes = 5 components — 3 gauge , 

2 tensor decaying modes + 3 velocity = 5 components . 



(A7) 
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